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Abstract 


We summarize our new results on cosmic strings presented in [1]. These results include: 
(l) The application of non-equilibrium statistical mechanics to cosmic string evolution. (2) 
A simple “one scale” model for the long strings which has a great deal of predictive power. 
(3) Results from large scale numerical simulations. (4) A discussion of the observational 
consequences of our results. An upper bound on G\i of approximately 10 -7 emerges from 
the millisecond pulsar gravity wave bound. We discuss how numerical uncertainties affect 
this. Any changes which weaken the bound would probably also give the long strings the 
dominant role in producing observational consequences. 
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We have recently concluded a major investigation of the evolution of cosmic strings 
using both analytical and numerical techniques. Our results are presented in detail in 
[1], and expand on work begun in [2]. The purpose of this letter is to present our most 
important results in a concise form. 

We consider only the simplest, directional and local strings, formed as defects at a 
unification phase transition. These accurately follow the Nambu equations soon after they 
are formed. In the simplest scenario, where the strings accrete matter gravitationally 
the observable consequences depend only on p , , the mass per unit length of the string[3]. 
This makes the theory highly predictive. As we shall explain, from improved numerical 
calculations and millisecond pulsar timing measurements, p. is now very tightly constrained. 

Vachaspati and Vilenkin[4] showed that the initial conditions for cosmic strings consist 
mainly of long strings, in the form of random walks. We take their initial configuration as 
our starting point. 

Statistical mechanics calculations [5] show that the initial state for the strings is 
rapidly pushed far out of equilibrium by the expansion of the universe. As the density in 
strings goes down, phase space heavily favors the long string chopping up into loops. One 
thus expects string interactions to work in the direction of equilibration, acting to speed 
up the rate at which the long string density pi falls. 

This process is well described by a “one scale model” [6] in which the length scale £, 
defined by 


( 1 ) 


PL , _ 

m i 2 

is assumed to be the only relevant scale in the problem. The long strings are random walks 
of step size 0(£) and mean separation 0(£). Equilibration proceeds through the chopping 
of small loops off the long string network. 

On dimensional grounds, the rate for this must take the form 


P'l = — - pl (due to chopping) 


( 2 ) 


Here p'l = and we use units where the speed of light is unity. The ‘chopping efficiency’ 
c is essentially a geometrical constant giving the rate, in units of £, that a random walk of 
step length £ chops itself into loops. One expects c to be a fairly small fraction of unity. 
We have tested (2) in exact flat spacetime simulations and find it accurately obeyed, with 
c fss 0.1 [7]. 
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In an expanding universe there are two other important effects. Hubble expansion 
decreases the string density and also stretches the long string. It is convenient to define a 
new variable 


Rh 

7 = T 


(3) 


being the ratio of the Hubble radius to the scale £ on the string. The full equation for this 
is [1] 

7 = -^(c7 ! -(2Bh-3)7-1) (4) 


The one-scale model has only one free parameter, c, and predicts the evolution of the 
long string density through the entire history of the universe. In the matter and radiation 
eras, Rjj = 3/2 and 2 respectively, and (4) has the stable fixed points 7 m = 
for c << 1. Whatever the initial string density, 7 approaches the appropriate fixed point 
and the network is ‘scaling’ (£ oc Rn)- Unlike previous analytic approaches [6], [8], our 
discussion of the statistical mechanics fixes the sign of c and makes it clear that the scaling 
solution is inevitable[6],[9]. 

In (1] we made a significant advance in numerical technique over previous work [2], 
[10], through the use of new nonsingular evolution variables. Our new crossing detection 
scheme is also a significant improvement over our previous very crude one [2]. 

We have found that with c = .074, equation (4) gives excellent agreement with all our 
expanding universe simulations in the radiation, matter, and transition eras. The impor- 
tance of the stable fixed point is that the state of the string network at late times depends 
very little on the precise initial conditions. This fact makes the resulting predictions more 
“fundamental”, and less susceptible to adjustments in the details of the theory. The fixed 
point behavior is also of practical importance. In the universe, the cosmic strings have to 
evolve for an expansion factor of more than 10 20 before matter accretion begins, whereas 
any present simulation can evolve them for a factor of 10 at best. 

It is important for many predictions to understand the evolution of loops chopped 
off the long strings. A single isolated loop can only break up into a finite number of 
offspring[ll]. Interactions between loops are necessary for equilibration into infinitesimal 
loops. However, the expansion of the universe increases the separation between loops and 
redshifts their velocities to zero, so all loops are eventually isolated and a distribution of 
non-intersecting loops remains. 
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We define PNl{l)dl to be the energy density in loops with ‘length’ l to l + dl, where l 
is the rest mass of the loop divided by p. The production of loops is described by 


PN/(0 = jlfNi(l/0 


(5) 


where }ni is a dimensionless function, the ‘energy production function’. Following the one 
scale model, the assumption is that at any time and for any string density the production 
of loops looks the same when everything is scaled according to We then define 


T n = f x n f NI (x)dx 

Jo 


( 6 ) 


to represent the different moments of fNi(x). 

One point we stress in [1] is that different observational consequences depend on 
different .F n ’s. For example, the zeroth moment To is related to c as follows from (2) and 
(5). The reason they are not strictly equal is that there is, in addition to the long string, 
an ‘intermediate’ population of intersecting loops with density pi « .5 pi. This results in 
To = c(pL + pi)/pi > c. 

The final non-intersecting loop distribution is given in the radiation era by 


and in the matter era by 


PNl(l) = 8 

Rpi 

(7) 

K = l^,/* 

(8) 

= R y 

(9) 


(10) 


In [1] we suggest that the process of loop emission may be viewed as ‘thermal’ radiation 
from the long strings. Even though the system is far from equilibrium, as with radiation 
from a hot body the spectrum of loops may be ‘thermal’. Detailed balance considerations 
then indicate a simple form for the energy production function 


f N i(x) - Ax ie 




(ii) 


for x > 1. 
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This form is well fitted by our numerical results, right down to our resolution with 
A « .27,5 « 2.6. 

The agreement of our results with the ‘thermal’ prediction, is intriguing but certainly 
no more than suggestive at this stage. There are many aspects of it we do not understand. 
For example, detailed balance yields the equilibrium production of all loops, whereas our 
results are for non-intersecting loops. 

We see good evidence that above our resolution the energy production function (5) is 
independent of time, as the one scale model predicts, and the various moments of /at/(*) 
-remain fairly constant throughout the simulation. However, in [1], we argue from the fact 
that chopping off must be favoured over reconnection that /jvj(x) must have an (integrable) 
x “5 divergence at small x in the radiation era scaling solution. Our results are consistent 
with this: we observe a buildup of /n/(®) at small scales as the simulation proceeds and 
the resolution improves. This contributes to the small scale resolution dependence. 

Since the various parameters defined above are all related to different moments of 
/at/(*)> the uncertainty in fsi gives us a rough measure of the uncertainties in the pa- 
rameters. How large are these? From higher resolution runs in the radiation era, and from 
flat spacetime simulations we have recently performed, we estimate our value of c could be 
too low by as much as a factor of 2. Considering the dependence on fyi, one finds from 
the above equations that 7 r could be too high by a factor of 2, A r too large by a factor of 
2? fa 3, 7 m could be too high by 25 and A m by a factor of 25. Here we have made the 
simplistic assumption that all Tn s have the same uncertainty. In fact the higher moments 
depend less on small scales and should have smaller uncertainties. We caution that there 
could also be systematic uncertainties such as those we discuss below which would alter 
these estimates. 

Recent results from the simulations by Bennett and Bouchet [12] differ qualitatively 
from ours in some respects. Their chopping efficiency is a factor of two larger, but that is 
consistent with our uncertainties. However, they see much more of their loop production 
all the way down at their scale of resolution. We believe the differences arise from the fact 
that our algorithm introduces numerical diffusion on the scale of our resolution (<< £) 
whereas the method of Bennett and Bouchet prevents the scale of kinks falling below a 
fixed minimal scale, of order their resolution. It is still not clear to us which method better 
approximates the continuum limit. 

Now we turn to the observational consequences of cosmic strings. The observation 
most likely to rule out the simplest gravitational accretion scenario is the millisecond 
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pulsar timing constraint [13], [14]. String loops radiate energy at a rate E = —TGp 2 with 
r typically 40-70 for simple loop trajectories [15], [16], [17]. If gravity waves with periods of 
years today are redshifted back to early times, one finds the waves were emitted around 10 6 
seconds, well before equal matter and radiation density. For these periods the density in 
gravity waves today may be calculated very simply from energy conservation and scaling. 
If pg(uj)du) is the energy density in waves of frequency u> to u + dw and p ra d the density in 
background radiation (photons plus neutrinos) then we find [1] 

. v Pg {v) = 2 2 x i 0 - 3 A r (G/z 6 )*r7 0 * (12) 

Prad 

where A r is given from (7), (8). Gp 6 = Gp x 10® and r 50 = T/50. 

From (12) we calculate the density in gravity waves divided by the density in the 
present microwave background p 7 = p ro< //( 1 + .23 N„) where N v is the number of neu- 
trino species still relativistic when the gravity waves were emitted. We take N v — 3. 
Unpublished results from millisecond pulsar timing indicate that today u)p g {fjj)j p 7 < 
3.2 x 10“ 3 [18], for waves with periods of 6 years, in the form of a stochastic background 
such as strings would produce. There are also weaker bounds at shorter periods - the 
maximum sensitivity is to waves whose period equals the total observation time. This 
gives the bound 

A P (G/i«)*r~* <0.9 (13) 

The observational limit is conservative: a full analysis has yet to be completed[l8]. In our 
simulations, A r ~ 10, which rules out strings with Gp > 10~®. However as we emphasized 
in [1], there are still large uncertainties. With our uncertainty estimates given above, we 
still rule out Gp > 10 -7 . This value is the minimum necessary for significant gravitational 
accretion by strings to occur by today [1], 

There are a few small effects that could weaken the bound. Late time splitting of 
loops or annihilation of string near cusps could conceivably reduce (12) by a factor of two. 
Therefore a reduction of the observational bound by a factor of 3 is still probably necessary 
to rule out the simplest cosmic string scenario. This should be possible soon[18]. 

Bennett and Bouchet have recently suggested a way around the problem [12] - from 
the fact that loop production evolves towards smaller and smaller scales in their latest 
simulations. They suggest that in the true scaling solution the long string may lose most 
of its energy to loops well below the resolution of current simulations. To see how small 
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these would have to be, one can show from the above that in the scaling solution A r is 
fixed by energy conservation to be 

X r = Y< > < U ) 

where <> means the expectation value weighted with / j\ri(x ) defined above, and A r is 
defined in terms of loop energies rather than masses. The quantity j r is as defined in (l),(2), 
where we include intersecting loops in the long string density. Bennett and Bouchet[12], 
and Allen and Shellard [19] report in recent simulations that 7 2 « 50. Even using this 
value, (13) and (14) still imply that < (//£)* > < -07. Thus to avoid the bound for 
G/ig = 1 typical loops must be produced with lengths of order 5 x 10 -3 £. With redshifting 
and late time effects one would still have to argue that most of the energy in the long 
string goes into loops of length one or two percent of the scale £. 

If this is the case, then it is clear that the dominant role in large scale structure 
formation would have to be played by the long strings. Detailed simulations of the accretion 
of matter around strings in the matter era are now underway [20]. It is already clear that 
long strings play a major role. 

In some situations, the gravity wave bound does not apply. Strings may couple to 
other long range fields and emit electromagnetic or goldstone boson radiation. For example 
the global string scenario [3] with Gfi = 10 -8 is still viable, though the dynamics of the 
network is much harder to calculate. 

We have also considered the gravitational lensing produced by strings in the matter 
era, and conclude that the fraction F of sky lensed by strings out to a redshift Z would be 

F*l<r 4 £ 2 GMA m + 7 L) ( 15 ) 

where, from our simulations we find A m = 2.1, and 7 m 4 = 7 m /4 = 1.0, with errors esti- 
mated above. The typical angle between lensed objects would be 2.6 Gfit arc seconds[21]. 
With our results, there would be approximately one lensed galaxy per 2 square degrees of 
the sky. Most galaxies would be lensed by long loops, and since galaxy clustering is quite 
weak at high Z, the lensed galaxies would typically be separated by a degree or so on 
the sky. Thus the Cowie-Hu event [22], where four apparently lensed galaxies were found 
within a region 40 axe seconds square, cannot be claimed as a successful prediction of the 
cosmic string scenario. 
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To conclude, we have reported several analytical and numerical advances in the under- 
standing of cosmic string networks. These advances have emerged from the application of 
statistical mechanics to cosmic strings, from the development of the one scale model, and 
from the use of our new non-singular variables in computations. We have also discussed 
some of the important observable consequences. The best chance of ruling out (or confirm- 
ing!) the simplest gravitational accretion cosmic string scenario lies with the millisecond 
pulsar. The present constraint appears to threaten the survival of this scenario, but a 
better understanding of the string network on very small scales will be necessary before 
there is general agreement on this point. More complex cosmic string scenarios, however, 
are not as tightly constrained. 
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